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Abstract 



We show that in realistic constituent quark models based on a two-body 
I confinement interaction color states appear in the middle of the experimen- 

ts ■ tally known spectrum. To avoid this situation we implement a three-body 



confinement interaction, introduced on an algebraic basis, into a semirela- 



<^ . tivistic version of the Goldstone boson exchange constituent quark model and 



by solving the Faddeev equations we show that this interaction can increase 
the gap between singlet and color states, such as the latter can be ignored 



and the known baryons can be described as simple q systems. We analyze 
ft, ■ the effect of a A- and a Y-shape three-body interaction. 
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I. INTRODUCTION 

In this study we are concerned with constituent quark models where the confinement 
interaction contains the two-body color operator Fi ■ Fj [1]. This means to assume that 
the one-gluon exchange picture survives to all distances. There is no real justification for 
applying this rule to the confinement, except simplicity. This relatively easy model has been 
widely used in baryon spectroscopy together with a hyperfine interaction generated either 
by the one-gluon exchange [2] or by the Goldstone boson (pseudoscalar meson) exchange 
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[3]. The Fi ■ Fj model has often been revived as being a convenient framework for studying 
q 2 q 2 systems (tetraquarks), see e. g. Ref. [4], q A q systems (pentaquarks), see e. g. Ref. [5] 
or q 6 system (the nucleon-nucleon problem), see e.g. Ref. [6]. 

The present work is devoted to baryons, seen as g 3 systems. The starting point is the 
observation that the presence of the Fi ■ Fj color operator implies that the Hamiltonian, 
besides color singlets, also possesses color octet and decuplet states at low energies. These 
states are undesirable. Their existence is in conflict with the experiment and with conclusions 
drawn from lattice or continuous formulations of QCD in the confining phase. In practice, 
these states have tacitly been ignored, relying on the assumption that they do not exist. 
However, as shown below, colored states are found close to singlets in the spectrum of 
consituent quark models. Therefore they should either be shifted to higher energies, by a 
change in the dynamics of the quark model, so that they could be ignored, or otherwise 
transformed into colorless states by coupling to gluons as for example in [7-9], or to color 
octet sea qq pairs as for example in [10]. 

A solution in the first direction has recently been proposed by Dmitrasinovic [11]. Al- 
though OCD has a local exact SUc(3) symmetry, in QCD-inspired potential models this 
symmetry appears as global, by construction. The observation was that by considering only 
two-body confining forces in a Hamiltonian model implies that the global SUc(3) symmetry 
is implemented only in a restricted way, namely the Hamiltonian is expressed in terms of 
the quadratic (Casimir) invariant operator of SU(3) only. More generally, one could express 
any quark model Hamiltonian inspired by QCD in terms of every SU(3) invariant operator. 
Based on this argument, in Ref. [11] a three-quark confining potential that depends on the 
cubic invariant operator of SU(3) has been added to the usual two-body confinement. In 
a qualitative way, it was shown that such an interaction can increase the gap between the 
singlet and color states, provided its strength has a specific sign and range. Besides its 
implications on the spectrum of ordinary q 3 systems, its role in exotic q 2 q 2 hadrons has also 
been considered. 

The three-body confinement interaction of Ref. [11] has subsequently been analysed in 
Ref. [12] in the context of a harmonic oscillator confinement where the size of its effect on the 
gap between the color and singlet states in the spectrum of a q 3 system has been estimated. 
Also the q 6 system with relevance to the nucleon-nucleon problem has been discussed. 

Below we first show explicitly how much the singlet and color states of a realistic con- 
stituent quark model appear close to each other. Next, by introducing a three-body confine- 
ment force related to the cubic invariant of SU(3), we demonstrate that one can separate 
them largely enough so that the color states can be ignored. For this purpose, here we con- 
sider the semirelativistic constituent quark model of Ref. [13], into which we implement a 
three-body confinement interaction. By using the Faddeev approach of Ref. [14] we calculate 
the lowest singlet, octet and decuplet states of a three-quark system, with this additional 
three-body confinement force. 

II. BASIS STATES 

Here we present the nucleon sector, namely S=l/2 and 1=1/2. The lowest color octet 
and decuplet q 3 states compatible with the Pauli principle are displayed in Tables I and 
II respectively. For orientation, their excitation energy E N — E , estimated in a two-body 
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harmonic oscillator confinement, is shown in the last column. In the same harmonic oscillator 
confinement the first colorless negative parity state is 395 MeV above the ground state. 
Therefore one can see that the + color octet and decuplet states appear substantially below 
this negative parity state. In the following we shall consider only the lowest color states, 
namely the + and V octets and the + decuplet. It is necessary to displace all color states 
at energies where they can safely be neglected, if the nucleon resonances are to be interpreted 
as pure g 3 excitations. This can indeed be achieved by an additional three-body force, as 
shown below. 



III. THE HAMILTON! AN 

The realistic total Hamiltonian presently under discussion has the form 

H = E K + ^) 1/2 + v *> + ^ + y x > (3- 1 ) 

i 

where pi and rrii are the momentum and the constituent mass of the quark i, V2& the two- 
body confinement, V^b the three-body confinement and V x the hyperfine interaction. Taking 
V 3b = one recovers the chiral constituent quark model of Ref. [13] which is the simplest 
semirelativistic version (no tensor, no spin-orbit) of the model originally proposed in Ref. 
[3], where the hyperfine interaction is spin and flavor dependent. The origin of such an 
interaction is thought to be the pseudoscalar meson exchange between quarks. This model, 
either in a nonrelativistic or a semirelativistic form, reproduces well the spectra of light 
baryons and in particular the correct order of positive and negative parity states. It has 
been shown that the addition of a tensor interaction to the spin-spin interaction in the 
Goldstone boson exchange model has a minor effect and thus maintains the quality of the 
spectrum [15]. Therefore, in the present study where the emphasis is made on the confining 
interaction the tensor force can be ignored. Although we have chosen a particular model, 
we expect our conclusions to be relevant for any reasonable constituent quark model with 
three valence quarks. The detailed parametrization of V x can be found in Ref. [13] and we 
do not reproduce it here, the main emphasis being laid on the confinement potential. We 
use a V 2 b term of the form 

Vn^v'irJil + FtFf), (3.2) 

i<j 

where = ^A" is the color charge operator of the quark i, \ a are the Gell-Mann matrices 
and 

"'fa) =V + -ynj , (3.3) 

with 7 = 2.33 fm -2 and Vo = —416 MeV as in Ref. [13]. The color part of V2& is consistent 
with Refs. [11,12] but different from that of [13]. It stems from the stability condition of a 
qq pair but it depends on an arbitrary constant which is fixed to 7/3. Anyhow, we rescale 
this interaction such as to reproduce the spectrum obtained in [13] (see below). 
The three-body confinement interaction has the form [11] 
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Vsb — Vijk — VijkCijk , 



(3.4) 



where Vijk is the radial part and Cijk is the color operator 

C ijk = d abc F? F) Ft , (3.5) 

with d abc some real constants, symmetric under any permutation of indices and defined by 
the anticommutator of the Gell-Mann matrices as {A a , A 6 } = 2d abc A c . The operator (3.5) 
can be rewritten in terms of the two independent invariant operators of SU(3) as 

r - - r _ -r^ + — ] (has 

where is the quadratic (Casimir) and the cubic invariant. The expectation values of 
(3.6) can thus easily be obtained from the eigenvalue of and [16] for any irreducible 
representation (A/x) of SU(3). These eigenvalues are 

(C {2) ) = \{\ 2 + /i 2 + A/i + 3A + 3/i) (3.7) 

and 

(C (3) ) = 4-(A - fj)(2\ + /J + 3)(A + 2(1 + 3). (3.8) 

Our attempt to choose a form for V^k is related to the present knowledge of confinement. 
Lattice calculations are ambiguous about the static three-quark potentials in baryons. Both 
Y-shapes [17] and A shapes [18,19] are supported. The difficulty to distinguish between them 
is due to the fact that the difference between the two gives at most 15 % in the dominant 
area law for the baryonic potential (for a recent review see e. g. [20]). Note that these 
results concern color singlet states only while here we deal with both color singlet and octet 
states. Potentials between static color sources in higher representations of SU(3), including 
the octet one, do also exist [21]. They satisfy the so-called Casimir scaling which means 
that potentials between sources in different representations are proportional to each other 
such as their ratios is given by the respective ratios of the eigenvalues of the corresponding 
(quadratic) Casimir operators. On the other hand the scaling proportional to the number 
of fundamental flux tubes embedded into a higher representation coincides with the Casimir 
scaling. In this situation we deemed reasonable to consider the working Ansatz that V^k 
has either a triangular shape 

V ijk = 7 c [\fi - fj\ + \ fj -r k \ + \r k - r*|] , (3.9) 

or a Y-shape 

Vijk = 1 c [\fl - f \ + \ fj - f \ + \ f k ~ fh\] , (3.10) 

with r the point where the three flux tubes meet such as to satisfy the SUc(3) gauge 
invariance (see e. g. [22]). Taking c = 1 and 7 = \/2y/a with a strength tension \[<j ps 1 
GeV fm _1 in (3.9), one recovers the linear confinement potential used in realistic consituent 
quark models as for example the model of Ref. [13] where the form (3.9) is used in conjunction 
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with a two-body color operator Fi ■ Fj which brings the factor 1/2 in 7. Taking c = 1 and 

7 = \fa in (3.10) one recovers the confinement potential used for example in Refs. [22], [23] 

or [24] where no color operator is present in the confinement part, thus one deals with colour 

singlets only. Here we use the radial form (3.9) or (3.10) in Eq. (3.4), i.e. in conjunction 

with a three-body color operator. In our calculations the value of 7 appearing both in (3.9) 

and (3.10) is taken to be the same as in Eq. (3.3) and the parameter c thus represents 

the relative strength of the three-body versus the two-body force. By assuming the same 

triangular shape both in V26 and V3& analytic calculations can be carried for a while. In this 

3 2 

way in Ref. [11] it was found that c must be located in the interval — — < c < -. The 

upper limit ensures that the lowest color singlet is below the lowest color octet state. The 
lower limit was required by the stability condition of the nucleon, (V2& + V3&) > 0. In Ref. 
[12] some arbitrariness was noticed for the lower bound because this is related to the choice 
of the color operator in (3.2). However, as the numerical results of Sec. IV will reveal, the 
above range of c is entirely satisfactory for our discussion of the triangular shape, because 
the gap of between the colorless and color states will turn out to be large enough. In fact 
the common conclusion of [11] and [12] was that c must be negative in order to obtain an 
increase in the gap between the color octet and singlet states in a g 3 system. Therefore in 
the calculations below related to the triangular shape, we take 

-- < c < . (3.11) 
2 v 1 

The Hamiltonian (3.1) was solved by using the Faddeev approach of Ref. [14], adequate 
for confining potentials. The necessary expectation values for the two-body color operator 
appearing in (3.2) are given in Table III. The expectation values of the three-body color 
operator (3.6) are taken from [11] or [12]. These are 10/9, -5/36 and 1/9 for the singlet, 
octet and decuplet SU(3) states, respectively. 

A particular advantage of the Faddeev approach is that the incorporation of permutation 
symmetry is very easy. For identical quarks, as in this case, the three Faddeev components of 
the three-quark wave function have the same functional form in their own Jacobi coordinate 
systems. Therefore the three equations can be reduced to a single one. From the structure of 
this equation it follows that the correct symmetry of the wave function under the exchange of 
any two particles is automatically guaranteed if the correct symmetry is implemented in one 
of the three components. In our calculation a bipolar harmonic basis was used, combined 
with spin, isospin and color basis states. If we select the basis states such as 

(_y+,+i+c = _ x ; ( 3 _ 12 ) 

where /, s, i and c are the relative angular momentum, spin, isospin and color quantum 
numbers of a two-quark pair, the Pauli principle for the three-quark system is satisfied. We 
have (— ) c = 1 when [f] c = [2] and (— ) c = - 1 when [f] c = [11], where [f] c is a given 
partition in the color space. If we denote the total angular momentum by L and the total 
parity by P = (— )' +A , where A is the relative angular momentum of the third particle with 
respect to the pair, then the lowest color octets must have L p = + , \~ and the lowest 
decuplet L p = + . In Table IV we show the structure of the lowest color states with 
the corresponding quantum numbers of their components. This treatment of permutation 
symmetry as well as the whole numerical procedure were checked against the exact harmonic 
oscillator results of Ref. [12]. 
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IV. RESULTS 



As mentioned above, our purpose is to understand the implications of a three-body color 
confinement interaction in a realistic model. We look separately at its effect either due to 
the A-shape (3.9), or to the Y-shape (3.10). 



A. The A-shape. 

This shape is currently used in conjunction with a two-body • Fj color operator as an 
approximation to the Y-shape [22]. Here we use it both in V2& and V^b as discussed above. 
Then the contribution of the color part of V2& and V3& sum up together to x% (i=l,8,10) as 
in Eq. (14) of Ref. [12]. 

+ — c i=l (singlet) 



Xi 



3 9 
13 5 

8 1 

- + -c i=10 (decuplet) 

3 9 



i=8 (octet) (4- 1 ) 



In Fig. 1 we show the dependence of some eigenvalues of the Hamiltonian (3.1) as a 
function of — c for —1.4 < c < 0. The case c = is the model [13] for which we reproduced 
the ground state nucleon mass 772^=940 MeV, and the resonance masses N(1440)l/2 + = 
1459 MeV, N(1535)l/2" = 1522 MeV and N(1710)l/2+= 1783 MeV. For c = the lowest 
color octets + and 1~ acquire the masses 1536 MeV and 1758 MeV respectively and the + 
color decuplet has 2077 MeV. Thus in a constituent quark model with three valence quarks 
the color states are so low that they cannot be ignored. If the three-body interaction (3.9) 
is switched on, the singlet states go substantially down and the color states go slightly up 
(or down) and a clear gap between the lowest color singlet S=l/2, 1=1/2 states and the 
color g 3 states emerges. This gap increases when c is decreased. At c ~ —1.4 it aquires a 
substantial value of more than 1500 MeV. 

Note however that the singlet states part of the spectrum shrinks and one has to bring 
it back close to the experimental spectrum. This can be achieved by a readjustment of the 
model parametres defined in Ref. [13]. In particular, in order to enlarge the distance between 
color singlet states one has to increase the strength tension 7 in (3.3) which will amount to 
an increase of the gap between the color and colorless states as well. It is then inferred that 
the gap will be large enough so that color states can safely be neglected. 

Before ending this subsection a word should be mentioned about the influence of the 
choice of the additional constant 7/3 in (3.2) on the results. As shown by Dmitrasinovic 
[11] the presence of this positive constant is crucial for the stability of a qq pair. In Ref. [12] 
some arbitrariness was noticed regarding this constant. In the g 3 spectrum, a change in the 
above additional constant amounts to a shift of the whole spectrum and some change in the 
slope of the confinement. If one wants to restore the spectrum to its initial form one has 
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to change the arbitrary constant Vq and the quantity 7 in Eq. (3.3). For example, if the 
additional constant is increased, 7 must be decreased. As a consequence the gap between 
the color and singlet states decreases. But this can be compensated by further decreasing 
the lower limit of c below - 1.5, which is allowed in the way discussed in [12]. 



B. The Y-shape. 

This is a genuine three-body force, both in the coordinate and the color space as well. 
Then the color part contributions of V2& and V3& does not sum up analytically to Xi as i n 
Eq. (4.1), because the expectation value of the radial part of V 2 b is different from that of 
Vzi,. Hence the numerical calculation of the entire spectrum is slightly more complicated. 

In the Y-shape, the flux tubes meet at 120° in order to ensure the minimum energy. If 
one of the interior angles of the flux-tube configuration is greater than 120° the minimum 
energy condition cannot be satisfied and the corresponding flux-tube collapses to a point 
which means that r*j= tq for i — 1, 2 or 3. Thus the Y shape moves continuously into a 
two-legged flux-tube configuration, where the legs meet at an angle greater than 120°. The 
geometrical arguments are clearly given in Ref. [22] for example. 

The spectrum associated to the Y-shape (3.10) is presented in Fig. 2. In this case the 
lower limit on c imposed by the inequality (3.11) is no more valid but c must remain negative. 
So we varied the parameter c between zero and - 2. At c = the spectrum is the same as 
in Fig. 1. 

The trend is similar to the A-shape case, but the mass of every colorless state vary more 
slowly with the strength |c| of the three-body Y-shape force. If one would rescale the |c|-axis 
by a factor 1/2, the A and the Y-shape results would look closer to each other. 

To better understand the role of the three-body confinement force we also calculated 
the root-mean-square radii associated to the states of Fig. 2. These are displayed in Fig. 3. 
One can see that the quark core radii of the ground state nucleon and its lowest N(J=l/2) 
resonances increase slowly with |c|. This is obviously the effect of the decrease of the con- 
finement contribution through the addition of the three-body term. On the other hand, the 
radii of the color states remain practicaly unchanged, similar to their energies. 

V. CONCLUSIONS 

Through the example of the Goldstone boson exchange model [13] we have shown that the 
spectrum of a realistic constituent quark model with a pairwise color confinement operator 
and a linearly increasing behaviour can accomodate both singlet and color low-lying q 3 
states. The lowest octet and decuplet color states appear in the middle of the observed 
spectrum which means they cannot be ignored. In this situation there are two alternatives: 

1) To keep simplicity. Then, as shown here, the addition of a three-body confinement 
interaction can increase the gap between singlet and color states by 1. 5 GeV or more. 
Then the color states can be neglected in calculations and the presently observed baryons 
can be described as pure q 3 states. To restore the quality of the spectrum one has to refit 
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the parameters of the model, including the strength tension 7, which has to be increased in 
order to enlarge the spacing between colorless states. The increase of 7 will lead to a further 
increase of the gap between the colorless and color states. 

2) If however the color states are maintained low in the spectrum of a Hamiltonian 
with a pairwise confinement interaction only, one has to give up the simple g 3 picture of 
baryons. The color states could give rise to colorless hybrid baryons having singlet-singlet 
+ octet-octet color components either of type (qqq)(qq) or of type (qqq)g. In the context 
of a potential model based on a one-gluon exchange hyperfine interaction, applied to the 
study of tetraquarks i. e. qqqq systems, it has been shown [4,25,26] that there is a strong 
mixing between states formed of two color singlet and two color octet (qq) pairs which leads 
to a substantial lowering of the variational energy. This should also be the case for hybrid 
baryons where singlet-singlet (qqq)i(qq)i and octet-octet {qqq)s{qq)8 components could mix 
substantially. The present study suggests that such hybrids are to be expected in all partial 
waves if only two-body Fj • Fj are considered. This is in agreement with the findings of Ref. 
[10] based on a two-body confinement and a schematic flavor and spin dependent hyperfine 
interaction model. The subject deserves a separate investigation inasmuch as the hybrid 
baryons raise a new interest [27]. 
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TABLES 



TABLE I. Three-quark color octet states [21] c compatible with the Pauli principle. The first 
column gives the orbital angular momentum L and parity P, the second, third and fourth columns 
give the permutation symmetry in the orbital, spin-isospin ad spin-isospin-color spaces respectively 
and the last column gives an estimates of the excitation energy En — Eq = (N + 3)/zu;[x8] 1//2 — 3fta; 
of each state for a harmonic oscillator two-body confinement with N quanta, oj/c = 2 fvrT 1 and 
Xs = Xs/Xi witn X8 and %i from Eq. (4.1) at c = 0. 
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[f]o 
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[f]isc 


N 


En — Eq 












(MeV) 
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[3] 


[21] 


[I 3 ] 
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1" 


[21] 


[3] 


[21] 
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1" 


[21] 


[21] 


[21] 


1 


616 


1~ 


[21] 


[I 3 ] 


[21] 


1 


616 


1+ 


[I 3 ] 


[21] 


[3] 


2 


1066 


TABLE 


II. Same as 


Table I but 


for color 


decuplet [3]c 


states where 


En — Eq = 


(N + 3)Mxio] 1/2 


— 3fkv with Xio = 


Xio/Xi where xw and xi are 


from Eq. (4.1) 


at c = 0. 
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N 
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[I 3 ] 





314 


1" 


[21] 


[21] 


[21] 


1 


814 


1+ 


[I 3 ] 


[3] 


[3] 


2 


1312 



TABLE III. The expectation values of the two-body color operator Oij = - + Fi ■ Fj of Eq. 

(3.2) (i < j) between three-quark states \{c2Cz)[f]c; ci[f]c) required in the Faddeev calculations. 
The particles 2 and 3 are first coupled to a symmetric [2] or an antisymmetric [l 2 ] state and next 
to particle 1 to a total color symmetry [f]c- 



[f]c 


[f]c 


Color operator 


Expectation value 
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[3] 




8/3 
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[11] 


[21] 
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[11] 
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O13 
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[11] 


[I 3 ] 


On 


5/3 
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TABLE IV. The quantum numbers of the components of a totally antisymmetric state 
of color symmetry [f]c, total angular momentum L, total spin S=l/2, and total isospin 1=1/2, 
compatible with the Pauli principle, as used in the Faddeev calculations. Here I, s, i and [f]c are 
the relative angular momentum, spin, isospin and color quantum numbers of a two-quark pair and 
A is the relative angular momentum of the third particle with respect to the pair. 
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FIG. 1. S=l/2, 1=1/2 eigenvalues of the Hamiltonian (3.1) with the A-shape three-body inter- 
action (3.9) as a function of — c. The solid line is the ground state nucleon, the dash-dotted line the 
N(1440) resonance (Roper), the double-dot-dashed line the lowest negative parity state N(1535), 
the double-dash-dotted line the N(1710) resonance. The color octets + and 1~ are represented by 
long- and short-dashed lines respectively and the + color decuplet by a dotted line. 
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FIG. 3. Root mean square radii of states of Fig. 2. 
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